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ABSTRACT 

We simulate the formation and chemodynamical evolution of 128 elliptical galaxies 
using a GRAPE-SPH code that includes various physical processes that are associ- 
ated with the formation of stellar systems: radiative cooling, star formation, feedback 
from Type II and la supernovae and stellar winds, and chemical enrichment. We find 
that the star formation timescale controls when and where stars form in the contract- 
ing gas cloud, determines the effective radius at given mass, and is constrained by 
observation to be ten times longer than the local dynamical timescale. We succeed 
in reproducing the observed global scaling relations under our CDM-based scenario, 
e.g., the Faber- Jackson relation, the Kormendy relation, and the fundamental plane. 
An intrinsic scatter exists along the fundamental plane, and the origin of this scatter 
lies in differences in merging history. Galaxies that undergo major merger events tend 
to have larger effective radii and fainter surface brightnesses, which result in larger 
masses, smaller surface brightnesses, and larger mass-to-light ratios. We can also re- 
produce the observed colour-magnitude and mass-metallicity relations, although the 
scatter is larger than observed. The scatter arises because feedback is not very ef- 
fective and star formation does not terminate completely in our simulations. ~ 25% 
of accreted baryons are blown away in the simulations, independent of the assumed 
star formation timescale and initial mass function. Most heavy elements end up locked 
into stars in the galaxy. The ejected metal fraction depends only on the star formation 
timescale, and is ~ 2% even to rapid star formation. 

Key words: methods: N-body simulations — galaxies: abundances — galaxies: el- 
liptical and lenticular, cD — galaxies: evolution — galaxies: formation 



1 INTRODUCTION 

The internal structure of galaxies, the spectrophotomet- 
ric, chemical, and dynamical properties at various locations 
within a galaxy, is determined by the processes of galaxy for- 
mation and evolution. Stars in a galaxy are fossils; the star 
formation and chemical enrichment history of the galaxy 
are imprinted on their kinematics and chemical abundances. 
The SAURON pro j ect with the William Herschel Telescope 
jBacon et alJl200ll:lEmsellem et alJl2004ft is providing wide- 
field mapping of the kinematics and stellar populations of 
nearby galaxies, which will certainly give stringent con- 
straints on galaxy formation and evolution. Multiobject and 
integral field spectrographs are being developed also on 8- 
10m ground-based telescopes, which will provide the time 
evolution of such internal structure. To infer the physical 
evolution processes from the observational data, it is neces- 
sary to construct a realistic model, i.e., a three-dimensional 
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chemodynamical model, and to compare the theoretical pre- 
dictions with such observational data. 

How elliptical galaxies form is a long-standing 
matter of debate. Two competing scenarios for 
the formation of elliptical galaxies ha ve bee n pro- 
posed; the monolithic collapse (e.g., iLarsonl 11974 



Arimoto & Yoshii 
Toomrd Il977t 



and the major merge r 



Kauffmann^WJiite^^^Md^rdon h993; 
Cole fc Frenkl Il99& ISteinmetz fc Navarro! 12002ft . 

hereafter K04), we constructed a 
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In Kobayashi (2004 
self-consistent three-dimensional chemodynamical model of 
ellipticals, introducing various physical processes associated 
with the formation of stellar systems; radiative cooling, star 
formation, feedback of Type II and la supernovae (SNe II 
and SNe la), and stellar winds (SWs), and chemical enrich- 
ment. We then argued that both formation processes should 
arise to explain the observed variation in radial metallicity 
gradients. The metallicity is enhanced in the central 
dense region, and the metallicity gradients are generated. 
However, because merging events weaken the metallicity 
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gradients, and because the secondary star burst induced by 
the mergers is not enough to regenerate them, galaxies that 
form monolithically have steeper gradients, while galaxies 
that undergo major mergers have shallower gradients. 
Therefore no correlation is found between the mass and 
gradi ents, as in the observation llKobavashi fc Arimotol 

While the internal structure of elliptical galaxies is 
greatly affected by their merging histories, their global prop- 
erties should be determined from their masses according to 
the scaling relations. Differences in merging history may pro- 
vide the scatter in these relations. In this paper, we inves- 
tigate whether our simulated galaxies follow the observed 
correlations; the Faber- Jackson relation, the Kormendy re- 
lation, the colour-magnitude relation, the mass-metallicity 
relation and the fundamental plane (FP). 

The FP is a c orrelation of early- type gal a xies with 2 + n 
param eters (e.g., iDiorgovski fc Davia Il987| iDressler et all 
1987) that reflect the internal structures, and is a clue to 
understand formation and evolution of early-type galaxies. 
One possible interpretation of the FP, defined by central 
velocity dispersions <7o, absolute effective radii r c , and sur- 
face brightnesses within an effective radius SB C attributes it 
to a correlation of the mass-to- light ratio Mj L to the total 
luminosity, or equ ivalently, to the total galaxy mass (e.g., 
iFaber et al1ll987t) . The dependence of M/L on the mass 
and luminosity stems from the stellar met allicity and/or 
age iPahre. Diorgovski fc de Carvalholll998l) . However, el- 
liptical galaxies may not be homologous along the FP, 
and there remains a dispersion that is not due to obser- 
vational error. This may be caused by metallicity, age, 
and/or dynamical disturbance. The F P is observed up to 
z ~ 0.5 in clusters <|Kelson et al.lll997L 2000), which is un- 
derstood as the evidence for passive evolution since z < 
1. For field galaxies, a comparable correlation is ob- 
served up to z ~ 1 llTreu et alJl200 it iGebhardt et alJl2003t 
l yande Ven. van Dokkum fc Franxll2003t Ivan der Wei et all 
2004). The larger scatter than in clusters and the zero-point 
offset at higher redshifts can be interpreted as an age differ- 
ence, although other possibilities such as dynamical distur- 
bance have not yet been discussed. 

The details of our GRAPE-SPH chemodynamical code 
were described in K04, and we briefly summarize them in 
321 Following the discussion in K04, we vary the parameters 
controlling star formation and the initial mass function, and 
we show how these affect global properties and their cor- 
relations f il3.2l . With the best parameter set, we show the 
scaling relations of simulated galaxies comparing with the 
observations. We focus on the fundamental plane, and dis- 
cuss the origin of the scatter in H3.3I §4 and §5 respectively 
contain the discussion and our conclusions. 



2 CHEMODYNAMICAL MODEL 

The characteristics of our GRAPE-SPH code may be sum- 
marized as follows (see K04 f or the detail ). 

i) The SPH method jMonaghanl Il992l for a re- 
view) is adopted, and the gravity is calculated in direct 
sum mation usin g the special purpose computer GRAPE 
iSugimoto et alJ Il990t . The SPH for mulation used in 
the code is almost the same as in iNavarro fc White! 



Table 1 . Number of simulated galaxies 







run 


ellipticals 


dwarfs 


cDs 


high resolution (TV 


~ 60000) 


13 


14 


4 





low resolution (TV < 


- 10000) 


52 


55 


37 





wider region (TV ~ 


60000) 


9 


5 


4 


9 


total 




74 


74 


15 


9 



Jl993l). The GRAPE-SPH code was originally written by 
Nak asato fc Nomotol i2003l) , and is highly adaptive in space 
and time through individual smoothing lengths and individ- 
ual timesteps. The calculations were done with the GRAPE5 
system in the National Astronomical Observatory of Japan 
and the GRAPE6 of the University of Tokyo. 

ii) Radiative cooling is computed using a metallicity- 
dependent cooling function. For primordial gas ([Fe/H] 
< —5), we compute the cooling rates using the two- 
bo dv_rjrocesses_of_jI_anri_He, and f ree-free emission, as 
in lKatz. Weinberg fc Hernauisli 1^996). For metal enriched 
gas ([Fe/H] ^ — 5), we use a metallicity-dependent cool- 
ing function computed with the MAPPINGS III software 
^Sutherland fc Dopital 119931) . In this cooling function, the 
elemental abundance ratios are set to be constant for given 
[Fe/H] according to the relations found in the solar neighbor- 
hood. [0/Fe]= 0.5 for Galactic halo stars for [Fe/H] < -1, 
and solar values for [Fe/H] ^ 0. We interpolate between 
these values for -1 < [Fe/H] < 0. 

i i i) Ou r star formation criteria are the same as in 
iKatzl lll992l) : (1) converging flow; (V-w)i < 0, (2) rapid cool- 
ing; t cool < t dyn , and (3) Jeans unstable gas; t dyn < t sound . 
The star formation timescale is proportional to the dynam- 
ical timescale (i s f = -tdyn), where c is our star formation 
times cale param eter. We also adopt the probability crite- 
rion dKatzlll992l) : A random number between and 1 is 
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compared with the probability P = 1 — exp 

time interval At s f = 2 Myr. 

If a gas particle satisfies the above star formation cri- 
teria, a fractional part of the mass of the gas particle turns 
into a star particle. Since an individual star particle has a 
mass of 10 5_7 Mq, it dose not represent a single star, but 
an association of many stars. The mass of the stars associ- 
ated with each star particle is distributed according to an 
initial mass function (IMF). We adopt a power-law IMF, 
4>{m) oc m~ x (the slope x — 1.35 gives the Salpeter IMF), 
which is invariant to time and metallicity. 

In H3.2I we discuss the dependence of our results on our 
free parameters; the star formation timescale c and the slope 
of the IMF x. In order to reproduce the observed radius- 
magnitude relation and mass-metallicity relation, we will 
choose c = 0.1 and x = 1.35 as a standard model. 

iv) For the feedback of energy and heavy elements, 
we do not adopt the instantaneous recycling approximation. 
Via SWs, SNe II, and SNe la, thermal energy and heavy ele- 
ments are ejected from an evolved star particle as functions 
of time, and are distributed to all surrounding gas particles 
out to a constant radius of 1 kpc. The ejected energy of each 
SW, SN H, and SN la are ~ 0.2 x 10 51 erg depending on 
metallicity, 1.4 x 10 51 erg, and 1.3 x 10 51 erg. We distribute 
this feedback energy in purely thermal form, although a frac- 
tion of it (given by a free parameter /kin) can be distributed 
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Table 2. Number of each types of galaxies with different star formation parameter c and the slope of the IMF x. 
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[El] 


[E2] 


[E3] 


[E4] 


[E5] 


[Dl] 


[D2] 


[D3] 


[D4] 


A: 


c=1.0, 


x=1.10 


5 


18 


19 


25 


11 


20 


13 


9 


4 


B: 


c=0.1, 


x=1.10 


10 


17 


15 


33 


16 


20 


18 


13 


2 


C: 


c=0.1, 


x=1.35 


8 


15 


12 


31 


20 


15 


1(5 


10 


2 



in k inetic form as a velocity perturbation to the gas particles 
(see lNavarro fc Whitdfl993f) . As shown in Fig. 14 of K04, if 
we adopt /kin = 0.1, the star formation efficiency is lower, 
the surface brightness of the final galaxy decreases at the 
centre, and metal-rich gas blows out, resulting in effective 
radii which are too large and metallicity gradients which are 
too shallow. 

For the metals, the mass-dependent nucleosynthesis 
yields of SNe II and SNe la are taken from Nomoto et al. 
1997ab. The upper and lower mass limits of the IMF are 
0.05 and 120M© , respectively. The progenitor mass ranges of 
SWs and SNe II are 8- 120M Q and 8- 5OM , respectively. 
For SNe la, we adopt the single degenerate scenario with 
the metallicity effect l)Kobavashi et al1ll998L 2000), where 
the progenitors are the Chandrasekhar WDs with an initial 
mass of 3 — 8Mq, and the lifetimes are determined from 
the lifetimes of the secondary stars with 0.9 — 1.5A/ Q and 
1.8 — 2.6Mq for the red-giants and main-sequence systems, 
respectively. 

v) The photometric evolution of a star particle is 
identical to the evolution of a simple stellar populat ion. 
Spectra f\ are taken from iKodama fc Arimotol (Il997t) as 
a function of age t and metallicity Z. 

vi) The initial condition is a slowly rotating sphere 
with a C DM initial fluctuati on generated by the COSMICS 
package feertschingerlll995l) . The cosmological parameters 
are set to be Ho = 50 km s _1 Mpc -1 , Q, m = 1.0, Q.a = 0, 
and as = 1.0. The initial angular momentum is added as 
rigid rotation with the constant spin parameter A ~ 0.02. 
We set different resolutions and total mass; the total mass 
of ~ 1O 12 M0 (baryon fraction of 0.1) with comoving radius 
of ~ 1.5 Mpc, and with N ~ 10000 and 60000 (the half for 
gas and the rest for dark matter), which are the same as in 
K04. The mass of a dark matter particle is ~ 1.8 x 10 8 M Q 
and ~ 3.0 x 1O 7 A#0, and the mass of a gas particle is ~ 

2.0 x 10 7 M Q and ~ 3.3 x 10 6 A/ Q , respectively. In this paper, 
we add a new sample for cD galaxies with the total mass of 
~ 10 13 Mq from wider initial conditions with radius ~ 3 
Mpc and N ~ 60000. This mass resolution is similar to the 
lower resolution of K04 sample. 

3 RESULTS 

3.1 EVOLUTION HISTORIES 

Simulating the chemodynamical evolution of 74 fields with 
different cosmological initial conditions, we obtain 128 
galaxies at the present time (i.e., t — 13 Gyr). Depending 
on the CDM initial fluctuations, in some cases one galaxy 
forms in the centre of the field, in the others one galaxy and 
several subgalaxies. We select galaxies having stellar masses 
in a 20 kpc sphere larger than 4.5 x 1O 7 M0. Although many 
less-massive subgalaxies form, we discard them because our 



resolution is not enough to study them in detail. We summa- 
rize the number of runs and the resulting galaxies in Table 
1. 

Different galaxies undergo different evolution histories. 
The difference is seeded in the initial conditions. Galaxies 
form through the successive merging of subgalaxies with 
various masses. These vary between a major merger at one 
extreme and a monolithic collapse of a slowly rotating gas 
cloud at the other. We classify galaxies into the 5 classes ac- 
cording to their merging histories, defining a major merger 
with / > 0.2 at z < 3; [El] monolithic, [E2] assembly, [E3] 
minor merger, [E4] major merger, and [E5] multiple major 
mergers (see K04 for the detail) . The numbers of galaxies of 
each class is slightly different from the adopted parameter 
set as summarized in Table 2, although the galaxy identifi- 
cation method is the same in K04. The fractions of major 
merger galaxies among elliptical galaxies are 46%, 54%, and 
59% respectively for the model A, B, and C. This is be- 
cause the number of merging events increase at 1 < z < 3 
with slower star formation (smaller c) and weaker feedback 
(smaller x). 

As well as the observed dwarfs, the variation of the star 
formation histories is seen for our simulated dwarfs with 
My, tot > — 19 mag. We classify these into the 4 classes; 
[Dl] initial starburst, [D2] continuous star formation, [D3] 
continuous star formation with recent starburst, and [D4] 
recent starburst. Observationally, the first class of galaxies 
would be dwarf ellipticals, the others are dwarf irregulars. 
The star formation rate (SFR) is truncated in [Dl] because 
of the supernova feedback. The intermittent SFR is induced 
by the gas accretion and/or interaction with other galaxies. 
In the following sections, we show the results both for giant 
and dwarf galaxies, but we should note that our numerical 
resolution in the dwarf galaxies is not enough; the particle 
number in a galaxy is small (~ 200 in the worst case) and 
the gravitational softening (1.0 kpc for the low resolution) 
is comparable to the size of galaxies. 

All simulated ellipticals form with an initial starburst 
&t z >2 with the typical timescale of 1 — 2 Gyr. The SFR de- 
creases because the gas is exhausted in the galaxy. The sec- 
ondary starburst is induced by the accretion of gas clumps 
and/or the merging of gas-rich galaxies. Not all merging 
events induce a secondary starburst; the fraction of merging 
events induce such a starburst is about 10%, depending on 
the gas mass of the secondary galaxy. The initial starburst 
is always larger than the secondary one. 

As discussed in K04, such a truncating SFR is due 
to the artificial cut-off of mass accretion caused by the 
vacuum boundary of our initial conditions. However, if we 
simulate wider region, the star formation continues longer, 
and thus colours tend to be too blue. The mass accretion 
can be continued but the star formation should have 
stopped by some process by z ~ 2 in cluster ellipticals (e.g., 
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Table 3. Total stellar mass, gas mass in the galaxy, gas mass outside the galaxy, the wind gas mass, and the total accreted baryon mass 
by the present in units of Mq . 





M«(< r 200 ) 


M eBS (< r- 2 oo) 


M gas (> 2r 20 o) 


M wind 


^acc 


-MwindA^acc 






A: 


4.8 x 10 10 


3.4 x 10 9 


4.1 x 10 10 


1.85 x 10 10 


7.4 x 10 10 


0.25 


0.39 


0.023 


B: 


5.6 x 10 10 


3.7 x 10 9 


3.6 x 10 10 


1.91 x 10 10 


8.1 x 10 10 


0.24 


0.34 


0.008 


C: 


5.9 x 10 10 


1.9 x 10 9 


3.5 x 10 10 


1.92 x 10 10 


8.3 x 10 10 


0.23 


0.33 


0.010 



Kodama & Arimotot Il997l lEllis et al.lll997l: iKodama et al. 


1998 


IStanford. Eisenhardt & DickinsorJll998j: iBrown et al. 


2onn| 


and z ~ 1 in field ellipticals (e.g., pilva & Bothun 


199c] 


ISchade et alJ Il999l: iBrinchmann & Ellis! l2000t 


Dade 


li. Cimatti & Renzinil I200C1: llm et al.l |2002fl. An anal- 



ogous truncation is required such as tidal stripping and 
effects of active galactic nuclei (AGN) (see Eventually, 
although the dark matter halo may be affected, the stellar 
population does not change so much. 

Contrary to the observational constraints, galactic 
winds are hard to generate in our simulations, and star for- 
mation does not completely stop by z = in most simulated 
elliptical galaxies. At the centre of the present-day galaxy, 
the dynamical potential is so deep that the gas density is 
high. In such regions, super metal-rich stars (Z ~ IOZq) 
keep on forming in the simulation. In several dwarf galaxies 
with M < 10 9 Mq, weak galactic winds can be seen even if 
/kin = 0. By the present epoch, roughly 25% of the baryons 
are blown away by the input of thermal energy from super- 
novae, but most of the heavy elements remain locked into 
stars in the galaxy (see H3.21 . 

3.2 PARAMETER DEPENDENCES 

We show the parameter dependence of the global proper- 
ties of a galaxy for varying star formation parameter c and 
slope of the IMF x in Figures and [5] The present val- 
ues are summarized in Table 3. We adopt identical initial 
condition for a non-major merger galaxy (ID 782389913, 
[El] monolithic-like) with high resolution. The model A 
(solid lines, c = 1.0, a; = 1.10) is the same as in K04, 
where the feedback is stronger than the others because the 
star formation takes place at earlier epoch before the gas 
cloud collapses enough. With the model B (dashed lines, 
c = 0.1, x = 1.10), the star formation takes place more 
slowly, and thus the stellar metallicity increases a lot. With 
the model C (dotted lines, c = 0.1, a; = 1.35), which is 
adopted in the following sections as a standard model, the 
stellar yield is reduced with smaller x, and thus the feed- 
back is weaker than in the other models. Here we define the 
galaxy as r < r2oo, i.e., the region with the mean density 
higher than 200p c at each time. 

Although the time evolution of the size and total mass 
(panels a and b) do not change at all, the stellar mass (panel 
c) varies with c and i by a factor of 20% and 5%, respec- 
tively, in the sense that stronger feedback decreases the stel- 
lar mass. The baryon fraction (panel d) is controlled by c. Al- 
though the dark matter contraction is independent of these 
parameters, c changes the peak epoch of the initial starburst 
(0.3 and 0.6 Gyr respectively for c = 1.0 and 0.1), which de- 
termines when and where stars form in the contracting gas 
cloud. On the other hand, the present gas fraction (panel e) 





Figure 1. Time evolution of global properties of a non-major 
merger galaxy; (a) the radius T200 with the mean density higher 
than 200p c , (b) the total mass M200 in r200i ( c ) the total stel- 
lar mass M t , (d) the baryon fraction b = (M» + M g )/Af200, (e) 
the gas fraction / g = M g /(M* + M g ), and (f) the mean stellar 
metallicity. The solid, dashed, and dotted-lines are respectively 
for models A, B, and C. 



does not depend on c, but on x. This is because the return 
fraction is smaller with x = 1.35 by a factor of 20 — 40% 
than for x = 1.10 depending on the lower mass limit Mi 
of the IMF; 0.33 (a; = 1.35, h'h = O.O5Af , model C), 0.42 
(1.35, O.1M ), 0.45 (1.10, O.O5A/ , model A), and 0.52 (1.10, 
O.IMq, model B). How much gas turns to stars finally de- 
pends both on c and x because some fraction of the gas is 
ejected from the galaxy. 

Although the residual gas mass in the galaxy mainly 
depends on x, its temperature is also affected by c. Figure [2] 
shows the time evolution of the gas mass in the galaxy and 
in the wind for all gas (panels a, d), hot gas with T > 10 6 
K (panels b, e), and cold gas with T < 10 4 K (panels c, f). 
With the strongest feedback (model A, solid lines), the half 
of the residual gas is hot, and no cold gas exists. With a 
longer timescale of star formation (model B, dashed lines), 
the residual gas can cool and the temperature decreases. 
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Figure 2. Time time evolution of the gas mass in the galaxy (a- 
c) and in the wind (d-f) for all gas (a, d), hot gas with T > 10 6 
K (b, e), and cold gas with T < 10 4 K (c, f). The wind gas is 
defined as the gas particles that have been inside r < r'200 an d arc 
outside r > 2r2oo a t present. The solid, dashed, and dotted-lines 
are respectively for models A, B, and C. 



With model C (dotted lines), the gas cooling is not efficient 
because of lower metallicity. 

In this simulation, although ten times larger amount 
of gas exists outside the galaxy (r > 2r2oo), some of them 
never accrete onto the galaxy. We define the wind gas as the 
gas particles that have been inside r < r^oo and are outside 
r > 2r2oo at present. The wind gas mass is almost the half 
of the outside gas mass. Although the temperature of the 
wind gas is different (panels e, f), the ejected wind mass 
by the present is independent of these parameters (panel d) 
with an ejected wind fraction, the ratio of the wind mass to 
the accreted baryon mass M w i n d/M acc , of ~ 25%. The wind 
efficiency, the ratio of the wind mass to the total stellar mass 
Mwind/M*, depends on c, and is ~ 0.4 and 0.35 for c = 1.0 
and 0.1. This means that it is easier for gas particles to be 
blown away when the gas collapse has not yet finished. 

The biggest effect of x is on the metallicity, and the 
mean stellar metallicity of the whole of the galaxy (panel 
f) is 2.2, 2.7, and Q.IZq respectively for model A, B, and 
C. T he observational estimate us ing the metallicity gradi- 
ents llKobavashi fc Arimotolll999l) is 0.5 — 1 solar. Although 
x = 1.35 is favored to meet this observational constraint 
under our star formation and feedback schemes, we cannot 
rule out other IMF. In principle, x = 1.10 gives stronger 
feedback, which can reduce the stellar metallicity. However, 
this feedback mechanism do not work well in our simulations 
and almost all the metals produced by supernovae are locked 
into stars. Because the metallicity-dependent cooling func- 




Figure 3. The star formation rates derived from the ages of stars 
in the present-day galaxy (top panel), the surface brightness pro- 
file (middle panel), and the oxygen abundance gradients (bottom 
panel). The solid, dashed, and dotted-lines are respectively for 
models A, B, and C. 



tion is included in our simulation, the metal enriched gas 
easily cools and loses the energy needed to escape. In spite 
of the strong dependence of metallicity on x, the ejected 
metal fraction (the ratio of metal mass in the wind to the 
total produced metals Mz, w ind/A/z,acc) depends only on c, 
and is 2% and 1% for c = 1.0 and 0.1, respectively. With 
larger c, metals start to be lost before the galaxy collapses 
enough, and thus more metals can be ejected. We should 
note, however, that these numbers are much smaller than 
expected from observations of X-ray gas in clusters. These 
suggest tha t the ejected m etal fraction should be larger than 
two-thirds iRenzinill2002Tl . 

As discussed in K04, the galaxies of the K04 simula- 
tion are too extended, which causes an offset in the radius- 
magnitude relation. This is because star formation takes 
place too early before the gas can accrete towards the centre. 
This problem can be solved by changing c. Figure [3] shows 
the star formation rates derived from the ages of stars in 
the present-day galaxy (top panel), the surface brightness 
profile (middle panel), and the oxygen abundance gradients 
(bottom panel). The peak epoch of the initial star burst is 
delayed from 0.3 Gyr with c = 1.0 (solid lines) to 0.6 Gyr 
with c = 0.1 (dashed and dotted lines). With c = 0.1, the 
surface brightness increases at the centre, which results in 
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A: c=1.0, x=1.10 B: c = 0.1, x=1.10 C: c = 0.1, x=1.35 




16 18 20 18 19 20 21 22 17 18 19 20 21 22 



- M v,tot[ ma g] -M vtot [mag] -M vtot [mag] 

Figure 4. The relations of effective radius r c (upper panels) and luminosity weighted metallicity (lower panels) to total V-magnitude for 
model A (left panels), B (middle panels), and C (right panels). The solid and dashed lines show the observed relations of giant ellipticals 
(Pahrc 1999) and dwarf ellipticals (Binggeli et al. 1984), respectively. See the caption of Fig[S]for the symbols, which show the merging 
histories for elliptical galaxies (larger points) and star formation histories for dwarf galaxies (smaller points). 



smaller r c (4.3, 2.7, and 2.2 kpc respectively for models A, 
B, and C). The metallicity gradient does not change. 



We show scaling relations for our simulated galaxies in 
Figure 3] The effective radius r c is derived by fitting the 
Vaucouleurs' law to the projection of \Z\ ^ 100 kpc on the 
X-Y plane (see K04 for the detail). Figure [I] shows the re- 
lations of effective radius r c (upper panels) and luminosity 
weighted metallicity (lower panels) to total V-magnitude. 
The solid and d ashed lines s how the observed relations for 
giant ellipticals (IPahrelll99Si) and dwarf ellipticals, respec- 
tively (Binggeli, Sandage & Tarenghi 1984). With model A 
(left panels), there is an off-set in the radius-magnitude re- 
lation because galaxies are too extended. The slope of the 
mass-metallicity relation is the same as observed, although 
the metallicity is too high. With model B (middle panels), 
the radius becomes smaller by a factor of three, which agrees 
well with the observed radius-magnitude relation. However, 
the metallicity is still too high and the slope of the mass- 
metallicity relation is too shallow. With model C (right pan- 
els), the radius-magnitude relation remains in good agree- 
ment with observations, and the metallicity decreases to 
meet the observations. Although the dispersion is larger, the 
simulated galaxies also follow the observed mass-metallicity 
relation. 



3.3 CORRELATIONS 

3.3.1 Scaling Relations 

For elliptical galaxies, it is well known that there are vari- 
ous correlations among physical scale parameters. The best 
known of these is the Faber- Jackson (1976) relation L tx <r n , 
where the slope is n ~ 4, but with a variation, depend ing 
on the sample definition jKormendv fc Diorgovskilll989t) . A 
correlation between the de Vaucouleurs' effective radius r e 
and surface brightness SB e was found bv lKormendvl il9T7f) . 
where more luminous galaxies have larger r c and fainter SB C . 
However, the scatter of SB C for a given a is quite large, and 
the SB C —a diagr am corresponds to a face-on view of the fun- 
damental plane feormendv fc Diorgovsk1ll989l) . For dwarf 
galaxies, these relations are different; the Mb — r c relation 
has a shallower slope, and the Mb — SB e re lation reverses 
with luminous dwarfs having brighter SB e iBinggeli et alJ 
Il984l) . 

Figure|S]shows the scaling relations between total stellar 
mass M, in r2oo, total V-band luminosity My, tot derived 
from the de Vaucouleurs' law, central velocity dispersion ctq 
in 2 kpc, effective radius r c , and mean surface brightness 
SB C within r c in the V-band. All points are for the simulated 
galaxies, and massive and luminous galaxies are on the right 
side in all diagrams. The symbols show the merging histories 
for elliptical galaxies and star formation histories for dwarf 
galaxies as listed in the figure caption. 

(panels a, b) The total stellar mass (r < 7-200) cor- 
relates well with the total luminosity derived from the de 
Vaucouleurs' fit (panel a) and with the central velocity dis- 
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Figure 5. The scaling relations among total stellar mass M* in r200, total V-band luminosity My, tot derived from the de Vaucouleurs' 
law, central velocity dispersion orj in 2 kpc, the effective radius r c , and mean surface brightness SB e within r c in V-band. All points are 
for the simulated galaxies, and massive and luminous galaxies are on the right side in all diagrams. The symbols, showing the merging 
histories for elliptical galaxies and star formation histories for dwarf galaxies, are: [El] monolithic (filled circles), [E2] assembly (filled 
squares), [E3] minor merger (filled triangles), [E4] major merger (open squares), and [E5] multiple major merger (open circles); [Dl] 
initial starburst (asterisks), [D2] continuous star formation (crosses), [D3] continuous star formation with recent star burst (plus), and 
[D4] recent starburst (three-pointed stars). The solid and dashed lines show the observed relations for early- type galaxies in the Coma 
cluster (Pahre 1999, also shown with the gray points) and for dwarf ellipticals in the Virgo cluster (Binggeli et al. 1984, B-V=0.9 is 
adopted) , respectively. 



persion (panel b) . The correlations for the stellar mass mea- 
sured in r e have a scatter which is twice as large. 

(panel c) The Faber-, Jacks on relation L oc a™ shows a 
smaller dispersion than observed (gray points) . The slope is 
n — 3.8 for this observation, but is steeper with n = 2.6 for 
the simulation. This is because the simulated dwarf galaxies 
tend to have smaller <tq because of lack of resolution. The 



slope for —21 < My < —19.5 mag is consistent with 2.7 
both for the simulation and the observation. 

(panels d, e, f) The mass-effective radius relation: Mas- 
sive/luminous ellipticals have larger effective radii. The sim- 
ulated galaxies follo w the observed relation in the Coma 
cluster jPahrd Il999l . solid line). The dispersion is not so 
small as ~ 0.7 dex, but is comparable to the observation 
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with ~ 0.5 dex (gray points). In the simulation, the surface 
brightness of the central part is smeared by the gravitational 
softening, which causes an uncertainty in fitting the de Vau- 
couleurs' law. For dwarf ellipticals, the observ ed relation has 
a sha llower slope than that of giant ellipticals feinggeli et alJ 
Il984l . dashed line). This tendency can be seen in the simu- 
lated dwarfs although the scatter is larger. 

(panels g, h, i) The mass-surface brightness relation: 
For giant ellipticals, massive/luminous galaxies tend to have 
smaller surface brightness, but this relation is reversed for 
dwarf ellipticals. These tendency can be reproduced in the 
simulation, but the dispersion is very large. 

(panel j) The surface brightness-effective radius rela- 
tion: The simulated giant ellipticals follow the observed 
relation where larger galaxies have lower surface bright- 
nesses in r c . The scatter is almost the same as observed. 
For dwarf galaxies, the observed relation is rectangular, and 
larger dwarfs have higher surface brightnesses. The simu- 
lated dwarf galaxies populate the same side of the observed 
relation, but the direction of the relation is different from 
observation. The effective radii of the simulated dwarfs tend 
to be too large, which is due to the lack of resolution. 

The observed scaling relations are reproduced in the 
simulations. The scatter exists even if the uncertainties of 
the simulations are taken account. The origin of the scatter 
is clearly demonstrated by the symbols. The galaxies that 
form monolithically (filled circles and squares) have smaller 
effective radii r G and thus brighter surface brightness SB C , 
while the galaxies that undergo major mergers (open circles 
and squares) have larger r c and thus fainter SB C . This is be- 
cause the merging events destroy the galaxy structures and 
make the radius larger and larger. The dynamical informa- 
tion is not fully wiped out but is blurred by a merger, as is 
shown using the difference of the energies of particles before 
and after merging events in Figs. 12 and 13 of K04. There- 
fore, we conclude that the scatter of the scaling relations 
stems from differences in merging histories. 

3.3.2 Mass-Metallicity Relation 

The colour-magnitude relation of elliptical galaxies is usually 
interpreted as luminous galaxies having higher stellar metal- 
licities. This is supported by the observation that a colour- 
magnitude relation with the sa me slope is found for high - 
redshift cluster ellipticals (e.g.. iKodama fc Arimotolll997t) . 
Previous studies showed that the line index (Mg2)o corre- 
late with the veloc ity dispersion ao at the galaxy centre (e.g., 
iDavies et alJl987l Bender, Burstein fc Faber 1993) and als o 
with total absolute magnitude (e.g., iBurstein et al]ll98St) . 
The relation with ao is tighter, but still shows some intrin- 
sic scatter. The same (Mg 2)o-oo relation is fou nd both for 
cluster and field ellipticals feernardi et al 11998ft . Therefore, 
the mass-metallicity relation of ellipticals is a common rela- 
tion independent of environment and time, which contains 
important information on the star formation and feedback 
processes during the early stages of galaxy formation. 

Figure HJ shows the sequences of stellar populations 
against the galaxy mass. The mass tracers are the total ab- 
solute V-magnitude and the central velocity dispersion. The 
characteristics of stellar populations are expressed by B-V 
and V-K colours, stellar abundances of oxygen [O/H] and 
iron [Fe/H], and stellar age, all of which are measured within 
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Figure 6. Colours of B-V and V-K (a-d), oxygen and iron abun- 
dances (e-h), and ages (i-j), which are measured in r c weighted by 
V-luminosities, versus total absolute V- magnitudes/central veloc- 
ity dispersions. The symbols are the same as in Fig[S] The solid 
line and gray points are for the observation taken from Bender et 
al. (1993) in (a-b) and Pahre (1999) in (c-f), respectively. 



r e and weighted by V-luminosities. Compared with observa- 
tion, the mass-metallicity relations (panels e-h) are weak, 
with shallower slope and larger scatter in the simulations, 
although the average is consistent. These are because the 
thermal feedback of supernovae is not enough to terminate 
star formation in the SPH simulations. Especially, ejection 
of the metal enriched gas do not work well. Although the 
steep slope of the observational relation requires not only 
more metal ejection in less-massive galaxies, but also more 
metal production in massive galaxies, such process cannot 
occur even if we change the IMF ( >I3.21 . For the scatter, 
we find no significant dependence on merging history; it is 
mainly caused by the age differences because the luminosity 
weighted metallicity is affected by bright young populations 
formed in the secondary starbursts. The relations for iron 
show larger scatter because iron is produced mainly by SNe 
la with longer lifetimes and affected by late star formation 
(sec [21 for the abundance ratios). 

Nonetheless, the majority of stars in the simulated gi- 
ant galaxies are formed in an initial starburst, and the 
luminosity-weighted ages are as old as 7 — 10 Gyr. No rela- 
tion is found between age and galaxy mass (panels i-j). For 
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dwarf galaxies, the ages decrease to 3 — 8 Gyr with a large 
scatter and a trend with mass can be seen, w hich is consis- 
tent with the observation using index l|Yamada et, all 
12004^1 . The scatter stems from the differing star formation 
histories of dwarfs. Among dwarf galaxies, dwarf ellipticals 
that formed by an initial star burst ([Dl], asterisks) have 
larger ages up to 5 — 8 Gyr. 

Therefore, the simulated colour-magnitude relations 
(panels a-d) have shallower slope and a larger scatter than 
the observations. At the reddest edge, we find the relation 
where massive ellipticals have red colours as observed. Red- 
der colours like V-K show smaller scatter because the origin 
of the scatter is the younger stars in the simulated galaxies. 

3.3.3 Mass-to-Light Ratios 

Figure |7| shows (a) baryon fractions, (b) gas fractions, (c) 
stellar mass-to-light ratios, and (d) total mass-to-light ratios 
against stellar masses. All of them are measured in spher- 
ical regions with the radius of 2r e , which includes ~ 60% 
of the total stellar mass. The baryon fractions (panel a) are 
~ 0.35 in massive galaxies, and decrease toward 0.1 in dwarf 
galaxies. In the galaxy, the baryon fraction increases toward 
the galaxy centre, but is 0.5 at the most in the central 2 
kpc. Even at the galaxy centre, equal amounts of dark mass 
exist with the baryons in the simulated galaxies. The gas 
fractions (panel b) are less than 5% for the galaxies with 
M» > 10 10 M Q , and increases towards 50% in dwarf galax- 
ies. Therefore, the stellar mass to the total mass ratio well 
correlate with the galaxy mass. More stars form in more 
massive systems. 

The stellar mass-to-light ratio (panel c) is derived in 
terms of the SSP model, which should be consistent with 
the observational estimates; for ellipticals, M*/L [Mq/Lq] 
~ 5 — 8 and 6 — 9 in the V-band and the B-band, respectively. 
A trend can be found that massive ellipticals have larger 
M t /L, although the scatter is large. For dwarfs, M*/L de- 
creases to 2 — 6 and 2 — 5 in the V-band and the B-band, 
respectively, which is due to younger ages and lower metal- 
licities than ellipticals. However, in our simulation, the con- 
tribution of dark matter is quite large even at the galaxy 
centre. The total mass-to-light ratios (panel d) are as large 
as 20 — 40 and 20 — 70 respectively for the simulated ellip- 
ticals and dwarfs. Observationally, the mass of dark matter 
can be estimated with the X-ray hot gas; the dark mat- 
ter mass can be several times larger than the stellar mass 
for several galaxies even at the galaxy ce ntre, which results 
in th e total mass-to-light ratios of ~ 20 iMatsushita et all 
1998). 



3.3.4 Fundamental Plane 

Figure El shows the V-band fundamental p lane shown in the 
K-space iBender. Burstein fc Fabeilll992f) . The parameters 
ki, K2, and «3 express masses, surface brightnesses, and 
mass-to-light ratios, respectively, and are defined as «i = 
(2 log 0-0 + logr e )/V2, k 2 = (21ogcr + 21og/ c - \ogr c )/y/E, 
and K3 = (2 log <7o — logic — logr c )/V3, where <to is the 
central velocity dispersion and I c = 10 
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solid line shows the observed relation for the V-band jPahrel 
Il999l) . and we reproduce the observed relations from the B- 
band to the near infrared. 
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Figure 7. (a) Baryon fractions, (b) gas fractions, (c) stellar mass- 
to-light ratios, and (d) total mass-to-light ratios against stellar 
masses, in spherical regions with the radius of 2r e . The symbols 
are the same as in FiglBl 



The K1-K2 diagram (lower panel) is the face-on view 
of the fundamental plane. There is no correlation between 
masses and surface brightnesses, and the simulated galaxies 
cover the similar regi on to the ob served giant and dwarf el- 
lipticals (gray points. IPahreHl999l ) . Dwarf galaxies populate 
the region with small masses and faint surface brightnesses. 

The K1-K3 diagram (upper panel) is the edge-on view. 
The simulated galaxies follow the observed re lation with a 
shallow slope (solid line, k 3 = 0.171/ti +0.143. IPahrJl99fll h 
where more massive ellipticals have large "mass-to-light ra- 
tios". Since the baryon fraction is as small as 0.5 even in 
the galaxy centre, the "mass-to-light ratios" expressed by 
K3 is affacted by the dark matter content. The rms fitting of 
the simulated ellipticals (dotted line), K3 = 0.247ki —0.035, 
is in agreement with the observation, although the slope is 
slightly steeper and the zero-point is larger by 0.1 dex than 
observed. Dwarf galaxies lie above the relation, where the 
total mass-to-light ratios are larger than for giant ellipticals 
(FigEJi). 

An intrinsic scatter exists along the fundamental plane. 
The origin of the scatter is clearly shown by the symbols; 
merger galaxies (open circles and squares) have smaller 
K2 and larger K3 than non-merger ellipticals (filled circles 
and squares). Figure [5] shows the histograms of the devia- 
tions from the fundamental plane for the non-major merger 
([E1]-[E3], gray area) and major merger galaxies ([E4]-[E5], 
hatched area). These distributions are different, and major 
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Figure 8. The fundamental plane shown in K-space for the V- 
band; the edge-on view (upper panel) and the face-on view (lower 
panel). The symbols are the same as in Fig|F] The solid line and 
gray points shows the observed relation (Pahre 1999). The dotted 
line show the rms fitting of the simulated ellipticals. 



merger galaxies have larger K3 than non-major merger galax- 
ies at given Ki. The thick dashed line is for the simulated 
dwarf galaxies, which have much larger K3. Therefore, the 
origin of the scatter along the fundamental plane is found 
to lie in differences in merging history. As discussed for Fig. 
|S] the galaxies that undergo major mergers tend to have 
larger r a and fainter I c . There is no significant change in 00 
and total luminosity L. From the definitions, these result in 
smaller K2 (oc L 2 aor~ 5 ) and larger K3 (oc L~ 1 a 2 r e ). 

What is the origin of the fundamental plane? In our 
simulation, the slope of the K1-K3 relation is originated from 
the combination of the metallicity, age, and the dark matter 
content, i) Metallicities: The slope mainly stems from the 
metallicity effect; metallicities are higher for massive galax- 
ies (FigE^f), which results in smaller total luminosities and 
larger stellar mass-to-light ratios (Fig[7J;) . The scatter of the 
mass-metallicity relation is not small, and the scatter of the 
FP can be reduced by the other effects, ii) Ages: No relation 
is found between age and mass (Fig|SJj), and galaxies are as 
old as 10 Gyr coevally at 3.5 < Ki <4. At m < 3.5, however, 
some low-mass metal-rich galaxies are younger and have K3 
as small as the other metal-poor galaxies (e.g., the galaxies 
at Ki — 3.25 and 3.13 on the dotted line), iii) Dark mat- 
ter content: Since the baryon fraction is larger for massive 
galaxies (FigJTJi), the total mass-to- light ratio is smaller for 
massive galaxies, which is the opposite of the FP. However, 
the old and metal-rich galaxies at K\ ~ 3.4 — 3.5 have larger 
baryon fraction and less dark matter content, which results 
in smaller K3 than the massive galaxies at ki ~ 4. 
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Figure 9. Histograms of deviations from the fundamental plane. 
The gray and hatched areas are for the simulated non-major 
merger ([E1]-[E3]) and major merger ([E4]-[E5]) galaxies, respec- 
tively. The thick dashed line is for the simulated dwarf galaxies. 



4 DISCUSSION 

The hydrodynamical simulation including star formation in- 
volves an uncertainty that is how to determine the star for- 
mation parameter c. We show here that it can be constrained 
from the scaling relation of galaxies. Because the star for- 
mation timescale controls when and where stars form in a 
contracting gas cloud in a dark matter halo, the star for- 
mation timescale determines not only the ages of stars, but 
also the size of the galaxy. To reproduce the observed radius- 
magnitude relation of elliptical galaxies, c is constrained to 
be 0.1 in our model, i.e., the local star formation timescale is 
ten times longer than the dynamical timescale. The global 
star formation timescale is found to be 1 — 2 Gyr, which 
is longer than the 0. 1 Gyr that is commonly adopted in 
one-zone models (e.g.. lKodama fc Arimotolll997t) . For spiral 
galaxies, the global star formation timescale should be sev- 
eral Gyr to meet various observations such as the metallicity 
distribution function of the Milky Way Galaxy. Physical pro- 
cesses that change the star formation timescale might be the 
existence of rotation, the suppression of star formation due 
to the UV background radiation, and environmental effects 
in clusters. 

We should note that mass accretion and star forma- 
tion are truncated artificially by the initial conditions in our 
simulations. However, in observed ellipticals, star formation 
should be truncated at z ~ 2 by some process, and an anal- 
ogous truncation is required, as discussed in H3.ll In our 
model, supernova feedback is not enough to stop the star 
formation, even if we increase the feedback energy with dif- 
ferent IMF ( iJ3.2H . Hypernovae, which eject ten times larger 
energy than normal supernovae, may increase the supernova 
feedback. Tidal stripping and ram-pressure stripping should 
affect, but it may be hard to explain the uniformity of cluster 
and field ellipticals. The AGN feedback can suppress the star 
formation effectively. If the rela tion between the black hole 
mass and the bulge luminosity iMaeorrian et al1ll99Sl) sug- 
gests that the AGN activity increases following the merging 
of galaxies, the AGN feedback may become effective around 
this redshift. 

To explain the lack of gas in present-day ellipticals, and 
to explain the heavy elements in the intracluster medium 
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JCiotti et alJll99ll) . a galactic wind seems indispensable. In 
our simulation, however, galactic winds do not occur in large 
galaxies, and thus star formation never terminates com- 
pletely. This causes the large scatter of 0.3 dex in the B-V 
colour- magnitude relation (Fig.|SJ. This problem arises from 
the SPH method and the feedback schem e. Including kinetic 
feedb ack /kin > as in other simulations llKawata fe Gibson] 
2003) does not seem to be a good solution. With /kin = 0.1, 
the surface brightness decreases at the centre of our ellip- 
ticals, and metal-rich gases blow out forming new stars at 
large radii. This results in effective radii which are too large 
and metallicity gradients which are too shallow to reproduce 
the observations. Changing the IMF slope is not good solu- 
tion either in our model, x = 1.10 gives stronger feedback, 
and the ejected wind mass does increase. Unfortunately, the 
ejected metal fraction does not increase, and the overall stel- 
lar metallicity becomes too high. The AGN feedback can not 
help solving problem because metals should be ejected dur- 
ing the star formation at very early epoch. It takes time to 
generate super massive black hole by the merging of galax- 
ies. Because of the metallicity dependence of gas cooling, 
enriched gas is easily turn to stars. 



Because of the existence of two distinct types of super- 
nova explosion that produce different elements on different 
timescales, the abundance ratios of the stellar population 
can be used to put constraints on star formation histories. 
SNe II, which are the core collapse-induced explosions of 
short-lived massive stars (> 8Mq), produce more O and 
Mg relative to Fe (i.e., [O/Fe] > 0) with a timescale of 10 6 " 8 
yr, while SNe la, which are the thermonuclear explosions of 
accreting white dwarfs in close binaries, produce mostly Fe 
and little O with a timescale of 0.5 — 20 Gyr. (The yields 
relative to solar value are almost the same between O and 
Mg both for SNe II and la.) Observationally, it has often 
bee n claimed that Mg is overabundant in elliptical g alax- 
ies jWorthev. Fa ber fe Gonzalez||l992|:|Thomas fe Marastonl 
120031) . Moreover. lFisher. Franx fe IllingworthHl995h showed 
that ellipticals with larger ao tend to have larger [Mg/Fe]o- 
However, in our simulation, luminosity weighted [O/Fe] 
spans —0.3 to —0.1 without any correlation with mass, al- 
though dwarf galaxies have larger [O/Fe] by 0.1 dex. (In 
model A with stronger feedback and the flatter IMF with 
x = —1.10, [O/Fe] spans —0.1 to 0.1.) Star formation needs 
to be terminated by efficient feedback in the simulations. 
We should note, however, that it may still be difficult to 
explain why [O/Fe] is larger in more massive galaxies. The 
timescale and duration of star formation should be longer 
for massive galaxies because the deep dynamical potential 
keeps gas cooling. So me possibilities have been suggested 
(e.g.. lMatteuccilll994l) ; i) the slope of the initial mass func- 
tion may be different in massive ellipticals, ii) the nucleosyn- 
thesis yields of SNe II may be different, namely, iron produc- 
tion may be different because of lower energies and/or larger 
fall-back, iii) the binary frequency may be smaller and less 
SNe la may occur in massive ellipticals, and iv) the metal 
enriched wind may cause selective mass loss, so that iron 
enriched gas can be ejected efficiently before it is consumed 
in forming the next generation of stars. 



5 CONCLUSIONS 

We study the formation and evolution of galaxies with a 
GRAPE-SPH chemodynamical model that includes various 
physical processes associated with the formation of stellar 
systems; radiative cooling, star formation, feedback from 
SNe II, SNe la, and SWs, and chemical enrichment. We 
simulate 74 slowly-rotating spherical fields with CDM ini- 
tial fluctuations (spin parameter A ~ 0.02), and obtain 128 
galaxies with stellar masses in the range 1O 9_12 M0 (74 ellip- 
ticals, 45 dwarfs, and 9 cD galaxies). In our scenario, galax- 
ies form through the successive merging of subgalaxies. The 
merging histories are various with differences seeded in the 
initial conditions. In some cases, galaxies form through the 
assembly of gas rich small galaxies, and the process looks 
like a monolithic collapse. In other cases, the final galaxies 
form through a major merger of preexisting galaxies. Ma- 
jor mergers are defined as those with mass ratio / > 0.2 
occurring at z < 3. 

Internal structure such as metallicity gradients is 
greatly affected by merging histories, while the global prop- 
erties are determined from overall masses according to the 
scaling relations. Assuming that the star formation timescale 
is ten times longer than the local dynamical timescale (i.e., 
c = 0.1), we succeed in reproducing the observed global scal- 
ing relations, e.g., the Faber- Jackson relation, the Kormendy 
relation, the colour-magnitude relation, the mass-metallicity 
relation and the fundamental plane. The different relations 
for ellipticals and dwarfs could be reproduced, although sim- 
ulated dwarfs have larger effective radii than observed be- 
cause of the lack of resolution. The luminosity- weighted ages 
of dwarfs span in wide range, 3 — 8 Gyr, depending on their 
star formation histories, while ellipticals are as old as 7 — 10 
Gyr independent of their mass. 

Adopting the Salpeter IMF (x — 1.35), we could re- 
produce the mass-metallicity relations both for the central 
stellar metallicity and for the mean stellar metallicity of the 
whole of the galaxy. However, the slope is shallower and 
the scatter is larger than observed, which are because the 
feedback is not so effective that most metals are locked into 
stars in the simulation. The colour-magnitude relation also 
shows a larger scatter because the star formation does not 
terminated completely in the simulations. 

An intrinsic scatter exists along the fundamental plane, 
and the origin of the scatter in the simulation lies in dif- 
ferences in merging history. Galaxies that undergo major 
mergers tend to have larger effective radii and fainter surface 
brightnesses, which result in larger ki (expressing masses), 
smaller K2 (surface brightnesses), and larger K3 (mass-to- 
light ratios). 

We examine the dependence of our results on the star 
formation parameter c and the slope x of the initial mass 
functions. Although the time evolution of the size and total 
mass do not change at all, the stellar mass depends on these 
parameters in the sense that stronger feedback (i.e., larger 
c and smaller x) decreases the stellar mass. We found that 
c controls when and where stars form in the contracting gas 
cloud, thus changing the baryon fraction and determining 
the effective radius at given mass. With the model to repro- 
duce the observed mass-radius relation, the baryon fractions 
are ~ 0.3 and ~ 0.1, and the total mass-to-light ratios are 
20 — 40 and 20 — 70Mq / Lq respectively for the simulated 
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ellipticals and dwarfs. The biggest effect of x is on the metal- 
licity, but it also changes the gas fraction and the residual 
gas mass in the galaxy. 

On the other hand, the wind gas mass, which is de- 
fined as the gas particles that have been inside r < r2oo 
and are outside r > 2r2oo at present, does not depend on 
these parameters, and ~ 25% of the accreted baryons can be 
blown away. The wind efficiency, the ratio of the wind gas 
mass to the total stellar mass, depends on c, and is 0.4 and 
0.35 for c = 1.0 and c = 0.1, respectively. However, most 
heavy elements end up locked into stars in the galaxy. The 
ejected metal fraction depends only on the star formation 
timescale, and is ~ 2% even if we take the quickest star for- 
mation rate (i.e., c = 1.0) under our scheme. To explain the 
metals detected in the intracluster medium in galaxy clus- 
ters, the feedback scheme should be improved so as that the 
enriched gas can blow away efficiently. Changing the IMF do 
not help in solving this problem under our star formation 
and feedback schemes. 
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